Astronomy & Astrophysics manuscript no. 2948 


February 2, 2008 


(DOI: will be inserted by hand later) 





The red optical afterglow of GRB 030725 

G. Pugliese 1 , P. M0ller\ J. Gorosabel 2 - 3 , B. L. Jensen 4 , J. P. U. Fynbo 4 , 
J. Hjorth 4 , S. F. J0rgensen 4 , B. Monard 5 , and C. Vinter 4 



in 
o 
o 



> 
(N 

in 
o 
in 
o 

Oh 
6 



1 European Southern Observatory (ESO), Karl-Schwarzschild-Str. 2, D-85748 Garching bei Munchen, Germany 

2 IAA-CSIC, PO Box 03004, 18080 Granada, Spain 

3 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

4 Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 K0benhavn 0, Denmark 

5 Bronberg Observatory, CBA Pretoria, PO Box 1 1426, Tiegerpoort 0056, Rep South Africa 

Received ; accepted 

Abstract. We present a photometric study of the optical counterpart of the long-duration Gamma Ray Burst (GRB) 030725, 
which triggered the HETE FREGATE and WXM instruments on July 25th, 2003, and lasted more than 160s. An optical 
counterpart was identified at the Bronberg Observatory in South Africa about 7 hours after the burst occurred. The optical 
afterglow (OA) was observed between 4 and 15 days after the burst with the 1.54m Danish telescope at La Silla in the V, R c , 
and I c bands. We fit a broken power law to the data and determine a break time in the light curve between 16 hours and 4.7 days 
after the first detection of the burst. The decay slope is a { = -0.59^' before and a 2 = -1.43 + 0.06 after the break. A bump 
may be present in the light curve, only significant at the 2cr level, 13.9 days after the main burst. The spectral slope of the OA, 
measured 12 days after the burst, is -2.9 ± 0.6 , i.e. it falls in the extreme red end of the distribution of previous OA spectral 
slopes. Observations of the field 8 months after the burst with the EMMI instrument on the NTT telescope (La Silla) resulted 
in an upper limit of R c =24.7 mag for the host galaxy of GRB 030725. The OA of GRB 030725 was discovered at a private, 
non-professional observatory and we point out that with the current suite of gamma ray satellites, an effort to organize future 
contributions of amateur observers may provide substantial help in GRB light curve follow up efforts. 
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1. Introduction 

Gamma Ray Bursts (GRBs) are powerful explosions in the 
y-ray band, with an emission which can peak up to a few 
MeV and a duration that varies between a few seconds up 
to hundreds of seconds (long GRBs). Since 1997, when 
the first optical counterpart was detected JCosta et al. 19971 
|van Paradijs et al. 1997) , several satellites have been launched 
to study GRBs and for half of the well localized GRBs, 
an optical afterglow (OA) has been discovered. As of to- 
day, there are about 60 detected GRB afterglows in the 
optical/IR bands. Detailed studies of GRB optical counter- 
parts provide unique information on the characteristics of the 
medium in which the afterglow evolves, the jet structure of 
the emitting region and the energy associated with the main 
burst. Nevertheless, a larger and well defined statistical sam- 
ple of well observed optical afterglows is necessary in or- 
der to address some still open questions concerning the uni- 
versal structure of the jet-like emitting region and a quasi- 
standard energy reservoir ( IFrail et al. 200U|Berger et al. 2 003 
IGhirlanda et al. 200 4 ), and the environment in which GRBs oc- 
cur (J^stetaOOOlll. 
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With the current suite of satellites, HETE-2, INTEGRAL 
and Swift, the rate of detected GRBs will reach an average of 
about 1 GRB every 3 days. Unless, however, the number of 
hours allocated to ground based optical follow-up will show a 
matching increase (about a factor of 10), the build-up of the sta- 
tistical sample will not proceed with the necessary speed, and 
many suitable afterglows will pass unobserved. Such a large 
increase is not realistic, and we must either explore alternative 
solutions or accept the loss of this opportunity to significantly 
increase our statistical sample. 

One possible alternative is to appeal to non-professional as- 
tronomers. The commonly used term "amateur astronomers" 
does not properly describe the high quality of the equipment 
which is often available at non-professional observatories, and 
historically amateurs have made important contributions in the 
fields of Solar System and variable star studies. The large num- 
ber of existent small private observatories could make signifi- 
cant contributions to GRB research in the fields of identifica- 
tion and high frequency early light curve sampling, the latter 
requiring that a suitable calibration of their data is available. In 
such a scheme a statistical sample of well documented GRB 
light curves could be created faster, and the professional GRB 
community would be free to be more selective in their follow- 
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up strategy. GRB 030725 is the first GRB optical afterglow to 
be discovered from a non-professional observatory, and there- 
fore serves as evidence that there is a valuable potential in the 
collaboration between amateurs and professionals in this field. 

GRB 030725, a bright long hard burst, was detected by the 
FREGATE and WXM instruments on board of the HETE-2 
satellite on July 25th, 2003 dShirasaki et al. 2003> . It had a flu- 
ence of ~ 2 x 10~ 5 ergcirr 2 in the 30-400 keV band, and a 
fluence of «2 x 10~ 5 erg cirr 2 in the 7-30 keV band. It showed a 
two peak shape, each peak with a FRED (Fast Rise Exponential 
Decay)-like structure. The first peak lasted about ~40 s, and af- 
ter about 160 s the second peak occurred with a duration of~ 10 
s. No detection of an afterglow in the X-ray or radio bands has 
been reported. 

In this article we present the results of the photometric ana- 
lysis of the optical afterglow of GRB 030725. The data were 
collected between 0.30 and 14.87 days after the main burst. We 
also observed the same field 8 months after the GRB occurred 
and obtained a limit on the magnitude of the GRB 030725 host 
galaxy in the R c -band. 



2. Observations and data reduction 

The optical afterglow of GRB 030725, located at R.A.= 
20 /l 33 m 59 i .47 and Dec. = -50°40'56".26 (J2000), with an er- 
ror of +0.1", was identified with the Bronberg 0.3m telescope 
approximately 7 hours after the trigger (Monard 2003), and a 
total of three exposures were obtained during the first night. 
Table 2 provides a log of those observations as well as of the 
later observations obtained with the 1.54m Danish telescope 
(1.54mD) on La Silla dVinter et al. 2003t . Observations at the 
1.54mD commenced 3.88 days after the burst during a period 
of adverse seeing conditions, and continued over the following 
1 1 nights. 

The Bronberg telescope is equipped with a direct CCD 
camera (SBIG ST-7E) with 765 x 510 pixel CCD, providing an 
18 arcmin E-W x 12.5 arcmin N-S field of view, with a pixel 
scale of 1 .52 arcsec per pixel. The camera operates in unfiltered 
mode only. 

The DFOSC on the 1.54mD is equipped with the 
MAT/EEV CCD 44-82 "Ringo" camera, a pixel scale of 0.39 
arcsec/pixel, and a corresponding field size of 13.7' x 13.7'. 
Because the Bronberg data were unfiltered we needed to de- 
fine an accurate calibration sequence with the 1 .54mD in order 
to combine the two data sets. We therefore first describe the 
reduction and calibration of the 1.54mD data, but will subse- 
quently return to the 0.3m data. 

The attempt to detect the host galaxy of GRB 030725 was 
performed with the NTT telescope at La Silla. Fig. ^ shows 
(left panel and right upper panel), the field of the GRB 030725 
observed with the 1 .54mD telescope on La Silla 3 days after the 
burst occurred. The right lower panel shows the same field ob- 
served with the 3.5m ESO-NTT telescope on La Silla, 8 months 
after the trigger. 
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Fig. 1. Left: The field of the GRB 030725 OA observed with 
the 1.54mD telescope on La Silla. The numbered circles indi- 
cate the standard stars used for the photometry. Their magni- 
tudes are shown in Table 1. Upper right: Zoom of the field 
(marked box in the left image) around the GRB 030725 OA 
observed with the 1.54mD telescope on La Silla 3 days after 
the burst occurred. OA marked with a circle. Lower right: The 
3.5m ESO-NTT image of the field of the afterglow obtained 8 
months after the trigger. Neither OA nor host galaxy is visible, 
its position is indicated with a circle. 

2.1. 1.54mD data reduction and calibration 

The data from the 1.54mD telescope in the V, R c , and I c bands 
were processed using standard IRAF 1 tasks for bias-subtraction 
and flat-fielding. The photometry was performed using the 
DAOPHOT package in IRAF (IStetson 1987IIStetson 19971 . 

We performed aperture and PSF photometry on both the 
optical afterglow (OA) and 7 isolated, bright unsaturated field 
stars which were used as reference stars for the relative pho- 
tometry. These stars are shown in the left panel of Fig.^ an d 
their corresponding magnitudes are shown in Table 1 . To cal- 
ibrate the light curve of the OA in the R c -band, two standard 
fields, MarkA and SA1 10, were chosen from the Landolt cata- 
logue of standard stars. These two fields were observed in the 
B, V, R c and I c bands. The corresponding calibrated magni- 
tudes of the optical afterglow in the V, R c and I c bands with the 
1.54mD telescope are shown in Table 2. 

2.2. The Bronberg data 

The first detection of the OT of GRB 030725 was obtained 
seven hours after the main burst occurred using unfiltered 
CCD imaging with an "amateur" telescope at the Bronberg 
Observatory. No photometric standards were obtained with the 
same instrument. Potentially the very wide response of unfil- 
tered data could introduce severe colour terms which would 

1 IRAF is distributed by the National Optical Astronomical 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. Magnitudes of the reference (ref) stars used for the photometry in the R c band. 



refl ref2 ref3 ref4 ref5 ref6 ref7 



Mag 16.12 


16.31 


17.97 


16.78 


17.96 


16.20 


19.08 


Mag.Error 0.04 


0.03 


0.04 


0.04 


0.03 


0.04 


0.03 



Table 2. Observing log (left) and photometric results (right) of the GRB030725 optical afterglow. "Unf" refers to the unfiltered 
data taken with the 0.3m telescope, while t - to refers to the observing epoch. The first 5 data points in the table were collected 
in July 2003, while all the following data were collected in August 2003. 



Date 


ExpTime 


Seeing 


Telescope 


t - 1 


Filter 


Magnitude 


Mag. Error 


(UTtime) 


(sec) 


(arcsec) 




(days) 








25.7944 


360 


3.1 


0.3m 


0.30 


Unf 


18.59' 


0.21 


25.8393 


360 


2.5 


0.3m 


0.38 


Unf 


18.67 1 


0.21 


25.9437 


360 


3.0 


0.3m 


0.48 


Unf 


18.90' 


0.21 


29.3355 


5500 


3.9 


1.54mD 


3.88 


Rc 


21.20 


0.04 


30.3151 


1800 


4.4 


1.54mD 


4.83 


Rc 


21.51 


0.05 


01.4347 


1200 


4.9 


1.54mD 


6.93 


Rc 


22.46 


0.55 


02.2920 


7800 


4.9 


1.54mD 


7.85 


Rc 


22.25 


0.06 


03.3268 


7800 


4.8 


1.54mD 


8.87 


Rc 


22.45 


0.08 


04.7899 


13800 


5.7 


1.54mD 


10.34 


Rc 


22.57 


0.10 


06.3535 


3600 


5.4 


1.54mD 


11.88 


Rc 


22.96 


0.10 


07.3176 


4200 


4.6 


1.54mD 


12.87 


Rc 


23.04 


0.09 


08.3116 


7800 


5.0 


1.54mD 


13.86 


Rc 


22.87 


0.17 


09.3360 


6600 


4.2 


1.54mD 


14.87 


Rc 


22.54 


0.20 


06.3211 


3600 


4.1 


1.54mD 


11.83 


V 


23.61 


0.15 


07.3383 


4800 


4.5 


1.54mD 


12.85 


F 


22.12 


0.17 


09.2839 


3600 


3.2 


1.54mD 


14.79 


F 


22.38 


0.19 



1 Unfiltered data calibrated to the R t .-band 



make it difficult to combine the unfiltered data with the later 
/? t -band data. However, the large gap in time between the 
Bronberg and 1 .54mD data renders the first three data points 
important in the determination of the decay parameters, and 
for this reason we went through a detailed calibration proce- 
dure aimed at determination of zero-point and colour terms of 
the unfiltered data. 

We followed a procedure very similar to that used for the 
photometric analysis of GRB 020813 by |Liet al. (2003)1 We 
selected 16 isolated stars which were unsaturated in both the 3 
unfiltered images and the 4 short exposure data in the B, V, I c 
and R c bands. For each of these 16 stars the calibrated colours 
B-V, V-R c , and R c -I c , as well as the unfiltered colours B-Unf, 
R c -Unf, I c -Unf, and V-Unf were computed. It turned out that 
the smallest gradients were found for R c -Unf, indicating that 
the effective central wavelength of Unf is close to that of R c . 
Because the unfiltered colour R c - Unf showed to be rather in- 
dependent of each of the calibrated colours for all the selected 
stars, we concluded that the R c -band was the one which best 
reproduced the unfiltered data. 

The OA R c -band magnitudes inferred from the calibrated 
unfiltered data are shown in the upper part of Table 2. 



2.3. The host galaxy 

The study of GRB host galaxies is an important step to de- 
termine GRB redshifts and to study the environment in which 
GRBs occur ( Jakobs son et al. 2 003 ). 

An attempt to detect the host galaxy of GRB 030725 was 
done on March 22, 2003, about 8 months after the main burst. 
The field in which the GRB OA was first detected was observed 
with the ESO NTT 3.5m telescope on La Silla. We used the red 
arm of the EMMI instrument in the R c -band and obtained 3 
exposures, each of 600 seconds, with a seeing of 0.6 arcsec. 
The data were processed using standard bias-subtraction and 
flat-field techniques and the photometry was computed using 
the same procedure described in the previous paragraph. There 
is no detection of the host galaxy and we compute a 2-sigma 
upper limit of 24.7 magnitudes in the R c -band for the GRB 
host galaxy. 

3. Light curve and spectral slope 

Many theoretical models have been developed to de- 
scribe the spectral and temporal evolution of GRB af- 
terglows JSari et al. 19981 Puglies e'et al. 1999) see also 
Zhang and Meszaros 2004 for a comprehensive review). Each 
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of these models assumes that the afterglow is created when 
relativistic particles are accelerated at the shock front between 
the flow generated by the GRB explosion and the interstellar 
medium (ISM), and emit synchrotron radiation visible in 
different bands. The distribution of the relativistic electrons 
is assumed to be a power-law, N(y e ) dy e oc y~ p dy e , where 
y e is the Lorentz factor. One of the main parameters which 
varies among all the theoretical models is the value assumed 
for p, the exponent of the electron distribution, whose range is 
generally chosen between 1 and 2.5, l |van Paradijs et al.~0 Q0 1. 
Correspondingly, each model will obtain different power-laws 
for the decay of the afterglow, which is well represented by the 
formula F v oc fv 13 , with a,/3 < 0. 

3.1. The light curve: Comparison with theoretical 
models 

The light curve of GRB 030725 in the R c -band is shown in 
Fig. |2] The data points fall in two groups, where the sec- 
ond group (1.54mD data) covers the longest span in time, 
and has the smallest error bars. Except for a possible upturn, 
or "bump" (see below), seen in the last few data points, this 
part of the curve is well fitted by a power-law with slope 
ar2 - — 1.43 ± 0.06. The slope of the three early data points 
alone is less well determined (a\ = -0.59 ± 0.67) but this can 
be further constrained by considering the break time. 

It is clear that a break must have occurred, but we are left 
with two possibilities. In case the break occurred earlier than 
3.88 days after the burst then the early power-law is the one 
stated above, with the additional limit that a\ > -0.98. If the 
break occurred later than 3.88 then the first 1.54mD data point 
falls within the pre-break light curve, and we find a\ -0.99 + 
0.04. Both scenarios are consistent with the data and both fits 
are shown in Fig. 3. We conclude that the break occurred before 
4.7 days, and that the early power law had a slope flatter than 
-1.03, i.e. a\ = — 0.59*2t?. For the purpose of determining the 
break time we only consider negative values of a\ . 

Breaks have been observed in the light curve of many after- 
glows (see Piran 2004 for an extensive review), and one of the 
possible explanations is related to an emitting region with a jet- 
like structure. According to this model, as long as 1 /F < 0j et , 
where 1 /T and 6j et are the relativistic beaming factor and the 
opening angle of the emitting region respectively, it will not 
be possible to distinguish between a spherical or a collimated 
emitting region. But when 1/T = 6j et , then the jet structure 
will become visible and this will correspond to a break in the 
power-law decay of the light curve. 

In the case of GRB 030725, the lack of data between 0.48 
and 3.88 days, did not allow an accurate estimate of the time 
of the break and the slope of the early power-law decay. Many 
models predict a change of slope between a few hours up to a 
few days after the main burst, consistent with our upper limit of 
4.7 days (observers frame). The allowed range of a\ reported 
above is in good agreement with slopes observed in other af- 
terglows, ac\ - -1, while the slope of the light curve after the 
break is in the extreme flat end of the distribution of other af- 
terglows -2.5 < a 2 < -1.5, J Andersen et al. 2000> . 




0.1 1.0 10.0 100.0 

(t-tO) days 



Fig. 2. Light curve of GRB 030725 in the R c -band. The solid 
line represents the fit of the data from the 1 .54 Danish telescope 
(diamonds) with a power law, assuming a break in the light 
curve 3.88+043 d avs a ^ ter tne burst and no SN contribution. The 
corresponding slope is equal to aq. -1.43 ± 0.06. The dashed 
line represents the fit of the first four data points with a power 
law and a corresponding slope equal to a\ -0.99 + 0.04. The 
dotted-dashed line represents the fit of the first 3 data points 
with a power law fitting with a slope of a\ - -0.59 + 0.67 and 
an earlier break in the light curve, 1 .30*?^ days after the burst. 
The arrow represents the upper limit to the magnitude of the 
host galaxy from the NTT data. 



If the break in the light curve of GRB 030725 is related 
to the presence of a collimated outflow, then the jet structure 
appeared between 16 hours and 4.7 days after the main burst, 
and p, the exponent of the power law distribution of electrons, 
should be less than 2. The lack of data during the first few days 
after the burst did not allow us to give a more accurate interpre- 
tation of this break. 

3.2. The light curve: A late time bump? 

Inspection of Fig. [5] shows suggestive evidence of an "up- 
turn" in the light curve after 13 days (observer frame). Such 
late time bumps in the light curves have been observed 
on several occasions, e.g. GRB 980326 which showed a 
red bump about 3 weeks after the burst (Bloo m et al. 19 99 1: 
JCastro-Tirado and Gorosabel 1999L GRB 011121 with a 
bump whose peak occurred between 13 and 75 days after the 
first y-ray observation dBloom et al. 2002> , GRB 020405 with 
evidence of a red late bump between 10 and 25 days after the 
main burst UPrice et al. 2003K GRB 02121 1 which showed the 
bump 27 days after the main burst ( Delia V alle et al. 20 03 ), and 
GRB 041006 with evidence of a bump 4 days after the burst 
( Sta nek et al. 20051 see lZeh et al. 2004l for a comprehensive re- 
view). For these GRBs the bump features were identified as SN 
contributions. 
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To test if our data support the presence of a SN we per- 
formed minimum x 1 two component model fits to the light 
curve where one component was a power law (the GRB) and 
the other was one of three historical SN light curves red- 
shifted to a se t of different redshifts (SN 19941, SN1998bw, 
and SN2002ap ( |Soderberg et al. 2005) ). None of the models in- 
cluding SN templates provided a significantly better fit than the 
fit of a pure power law, so we conclude that there is no signifi- 
cant evidence that the marginal (2 <x) upturn after 13 days can 
be identified with a SN bump. 

3.3. Spectral slope 

The OT was detected in three bands (V, R c , I c ) 12-13 days after 
the burst. On the assumption that the ct2 decay was monochro- 
matic, we computed (including corrections for Galactic red- 
dening), V-R c = 0.70+°^ and R c - I c = 0.87+° ^ (two sided 
errors computed using ML techniques, see Appendix A in 
|Fynbo et al. 20 02 1. In order to compute the spectral slope /? 
(F v oc \P, f} < 0) we transformed the V, R c , I c magnitudes 
to the AB system using the offsets to the AB system given in 
|Fukugita et al. (1995 ) and fit a linear relation to corresponding 
values of v c and AB. Here v c is the central frequency for the 
filters. In this way we infer a spectral slope of /3 — -2.9 + 0.6. 
Even allowing for the rather large error bar, this value is in the 
extreme red end of the distribution of spectral slopes, where 
the typical value is around 1 and a few limit cases around 2 
(Sim on et al. 200ll . We note that the OT colour was obtained 
12 days post-burst, so we can interpret it as either the result of 
significant extinction due to dust in the host galaxy, or the result 
of a contribution from an underlying SN, or the GRB occurred 
at a high redshift. 

4. Summary and discussion 

Since the discovery of the first GRB optical coun- 
terpart in 1997, about 300 GRBs have been detected 
by several satellites, but only for about 25% of them 
has an optical afterglow been detected (from webpage: 
http://www.mpe.mpg.de/~jcg/grbgen.htmli. For even smaller 
percentage, about 15% of the total detected GRBs, the redshift 
has been determined and only a few of these showed a late time 
bump in the optical light curve. 

GRB 030725 was a long hard GRB with a slope a x = 
-0.59^44 of the early time optical light curve, a break be- 
tween 16 hours and 4.7 days after the main burst, and a slope 
ai = -1.43 + 0.06 after the break. On the assumption that 
the break is related to a collimated outflow, and using the re- 
lation between the opening jet angle and the time in which the 
break in the light curve occurs ( |Pugliese et al. 1 999 1, we find 
jet opening angles of 7° and 12°, at 16 hours and 4.7 days re- 
spectively. 

We detected a weak feature of a late time upturn in the light 
curve 14 days after the burst and tried to model it with the 
contribution from SNs at different redshifts, using SN 19941, 
SN1998bw, and SN2002ap as templates. We found that none of 
these SNs provided a better fit to the data than a simple power 
law. 



We did an attempt to detect the host galaxy of GRB 030725 
about 8 months after burst, and we derived a 2-sigma upper 
limit of R c =24.7 mag for the GRB host galaxy. 

We determined the colours of GRB 030725 (V - R c = 
0.70^ and R c - I c = 0.87+° ^) corresponding to a spec- 
tral slope of (3 — -2.9 + 0.6. T his is significantly red- 
der than the majority of GRB OAs (Sim on et al. 200ll . The 
characteristics of the environment in which GRB occur are 
still an open question, ( Berger et al. 2002 ] |Fynbo et al."2 003 
ICastro Ceron et al. 20041 ITanvir et al. 2004t . but so far no 
GRB OA was found with such a steep spectral slope. 
This redder slope could be due to either high redshift, 
or contribution from an underlying SN, or dust absorp- 
tion in the host gal axy. Some authors dHjorth et al~2 003 ; 
IProchaska et al. 20041 |Vreeswijk et al. 2004} have shown that 
some GRB host galaxies have a rather low amount of dust and a 
rather low metallicity. The environment of GRB 030725 could 
be very dense along the line of sight to the GRB itself. Further 
observations in IR to mm wavelength bands would be interest- 
ing as they would allow us to better understand the nature of 
this spectral slope. 
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